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I report on a lattice computation of the energy of a system of two light 
quarks and two static antiquarks as a function of the separation of the 
static antiquarks. In terms of hadrons such a system corresponds to a pair 
of B mesons and its energy to the hadronic potential. I present selected 
results for different isospin, spin and parity combinations of the individual 
B mesons mainly focusing on those channels relevant to determine, whether 
two B mesons may form a bound tetraquark state. 

PACS numbers: 12.38. Gc, 13.75.Lb, 14.40.Nd. 

1. Introduction 

Lattice computations of the potential of a pair of static-light mesons (in 
the following also referred to as B mesons) are of interest, because they con- 
stitute first principles determinations of a hadronic force. Such potentials 
can e.g. be used as input for phenomenological calculations to determine, 
whether two B mesons may form a bound tetraquark state. 

In the literature interactions between static-light mesons have been stud- 
ied in the quenched approximation [U El [3l [H [5] and recently also with dy- 
namical quarks [6l [7] . Here I report on the status of an investigation with 
two flavors of dynamical Wilson twisted mass quarks. Forces are not only 
studied between the lightest static-light mesons (denoted by S), but for the 
first time also first excitations are taken into account (denoted by P-)- 
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2. Trial states and quantum numbers 

Quantum numbers of single static-light mesons as well as of pairs of 
static-light mesons (BB systems) have been discussed in detail in In 
the following I give a brief summary. 

2.1. Static-light mesons 

I consider static-light mesons made from a static antiquark Q and a 
light quark G {u , d} without non-trivial gluonic excitations. They can 
be labeled by the z-component of isospin 1^ = ±1/2, the z-component of 
the light quark spin jz = ±1/2 and parity P = ±. The V = — static-light 
meson (denoted by S", corresponding to B/B* in [8]) is the lightest static- 
light meson. It is lighter by around 400 MeV than its parity partner with 
V = + (denoted by P_, corresponding to Bq and B^). 

The corresponding static-light meson trial states are Q^^tp\Q) and 
Q^jil)\Vt) for S mesons and QVl^) and Qtj^sV'I^) for P_ mesons, respec- 
tively. 

For a more detailed discussion of static-light mesons I refer to [9, 10]- 

2.2. BB systems 

The aim of this work is to compute the potential of a pair of B mesons as 
a function of their separation R (without loss of generality I choose the axis 
of separation to be the z axis). To this end one has to compute the energies 
of eigenstates of the QCD Hamiltonian containing two static antiquarks 
0(ri) and Q(r2), ri = (0,0,-i?/2) and ra = (0,0,+i?/2), which define the 
positions of the two B mesons, and which will be surrounded by light quarks 
and gluons. 

These BB states are characterized by the following five quantum num- 
bers: isospin / e {0 , 1}, the z-component of isospin 1^ € { — 1, 0, +1}, 
the absolute value of the z-component of the light quark spin \jz\ G {0 , 1}, 
parity 7-" = ± and "x-parity" (reflection along the x-axis) Vx = ^■ 

I use BB trial states 

(Cr)^B(Qc(ri)V?^(ri))(Qc(r2)V'g^(r2))|n), (1) 

where the lower indices A, B and C denote spinor indices, C = 7072 is 
the charge conjugation matrix and F is a suitably chosen combination of 7 
matrices. Note that it is essential to couple the light degrees of freedom of 
both mesons in spinor space, because these degrees of freedom determine 
the quantum number Proceeding in a naive way by coupling light 

and static degrees of freedom in both B mesons separately will not result 
in a well defined angular momentum and, therefore, will mix different 
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BB sectors. To obtain 7 = 0, the flavors of the hght quarks have to be 
chosen according to = ud — du, while for 1 = 1 three possibiHties 

exist, ip^-^^ ■ip^'^^ s {uu , dd , ud + du}. For a list of BB trial states and their 
quantum numbers I refer to f6], Table 1. 

3. Lattice setup 

I use 24^ X 48 gauge held configurations generated by the European 
Twisted Mass Collaboration (ETMC). The fermion action is Nj = 2 Wil- 
son twisted mass |1H [T2] at maximal twist, where static-light mass differ- 
ences are automatically 0{a) improved [9j. The gauge action is tree- level 
Symanzik improved [TH]. I use gauge coupling /3 = 3.9 and light quark 
mass fiq = 0.0040 corresponding to a lattice spacing a = 0.079(3) fm and a 
pion mass mps = 340(13) MeV [H]. For details regarding these gauge field 
configurations I refer to \15\ [T6] . 

In twisted mass lattice QCD at finite lattice spacing SU(2) isospin is 
explicitely broken to U(l), i.e. Iz is still a quantum number, but I is not. 
Moreover, parity V has to be replaced by twisted mass parity 

p(tm)^ which is 

parity combined with light fiavor exchange. The consequence is that there 
are only half as many BB sectors in twisted mass lattice QCD as there 
are in QCD, i.e. QCD BB sectors are pairwise combined. Nevertheless, it 
is possible to unambiguously interpret states obtained from twisted mass 
correlation functions in terms of QCD quantum numbers. The method has 
successfully been applied in the context of static-light mesons [Hj and is 
explained in detail for kaons and D mesons in [18]. For a more elaborate 
discussion of twisted mass symmetries in the context of BB systems I refer 
to an upcoming publication [19]. 

4. Selected numerical results 

The potential of a pair of B mesons is extracted from the exponential 
falloff of correlation functions of trial states ([T|). As explained in subsec- 
tion 12.21 the trial states differ by the spin coupling of the light quarks via 
the 4x4 matrix F (16 possibilities) and by their light quark flavor (4 possi- 
bilities). Consequently, there are 16 x 4 = 64 different correlation functions. 
This number can also be understood from the point of view of individual B 
mesons: since each of the two B mesons has 8 possibilities regarding quan- 
tum numbers {Iz = ±1/2, jz = ±1/2, V = ±), there should be 8 x 8 = 64 
BB correlation functions. 

As outlined in subsection 12.21 the BB trial states ([1]) can be classified 
according to BB quantum numbers. However, to interpret the BB poten- 
tial obtained from the correlation function of a specific trial state ([T|), it is 
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also useful to express this trial state in terms of individual B mesons, e.g. 
^(1)^(2) =ud-du , T = j5 {-S-tSi + S^S-^ - P^Pi + PiPt)\^) 

V,(i)V;(2) =ud-du,T = 7075 ^ i-S^S^ + S^S^ + P^P^ - P^P^m), 

where S and P denote creation operators for S and P_ mesons, respec- 
tively, and t and J, indicate the orientation of the light quark spin. After 
linearly combining the two trial states via F = 75 + 7075 to eliminate the 
P_P_ contributions, one can extract an SS potential with quantum num- 
bers {I, Iz,\jz\,P,Px) = (0,0,0,—,+). Similarly, one can estimate a P^P^ 
potential with the same quantum numbers by choosing F = 75 — 7075. Re- 
sults are shown in Figure [T]^a). Further examples are discussed in [6]. 



(a) attractive SS and P_P_ singlet 
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(c) isospin triplet (1=1 , |j^|=0) 
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(b) singlet (1=0, ljj=0) 
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(d) isospin/z-rotatlon sextet (1=1, |j^|=1) 
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Fig. 1. BB potentials as functions of the meson separation, (a) SS potential and 
P^P^ potential both with quantum numbers {I,Iz,\jz\,P,Px) — (0,0,0,—,+). 
(b), (c), (d) Attractive SS singlet, isospin triplet and isospin/spin sextet. 



BB potentials mainly differ in their asymptotic value at large meson 
separations, which is approximately 2m{S), m(S) + m{P-) or 2m(P_) for 
SS, SP- or P-P^ combinations, respectively, and whether they are attrac- 
tive or repulsive at small meson separations. Some of the potentials, even 
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though they differ in their quantum numbers, are exactly degenerate due to 
isospin symmetry {1=1 triplets) or rotational symmetry around the ^;-axis 
(Ij'^I = 1 doublets). In summary the number of attractive and repulsive BB 
potentials for SS, SP- and P^P^ combinations and their degeneracies are 
as follows: 

SS potentials, attractive: 1 © 3 © 6 (10 states). 

repulsive: 1©3©2 ( 6 states). 

SP- potentials, attractive: 1©1©3©3©2©6 (16 states). 

repulsive: 1©1©3©3©2©6 (16 states). 

P-P- potentials: identical to SS potentials (16 states). 

Consequently, from the 64 trial states ([TJ one can extract 24 different po- 
tentials. 

Of particular interest, when asking, whether two B mesons may form 
a bound tetraquark state, are attractive SS potentials. As stated above 
there are three different possibilities, a singlet, an isospin triplet and an 
isospin/spin sextet. The three potentials are shown in Figure mb) to^d). 
While the triplet and the sextet are rather similar, the singlet is significantly 
different: it is deeper and wider and, therefore, probably the best candidate 
to start investigations, whether there are bound BB tetraquark states. Such 
phenomenological investigations will be part of a subsequent publication. 



5. Conclusions 

I have presented selected results of a first principles lattice computation 
of BB potentials. Various channels characterized by the quantum numbers 
{I,Iz.i\jz\-,'P-,'Px) have been investigated. Compared to existing publica- 
tions the computations have been performed with rather light dynamical 
quark mass (mps ~ 340 MeV) . The results have been interpreted in terms 
of individual S and P- mesons. Future plans include studying the light 
quark mass dependence, the continuum limit as well as finite volume ef- 
fects. Moreover, also BBg and BgBg potentials could be computed. To 
treat the s quark as a fully dynamical quark, such computations should be 
performed on Nf = 2 + 1 + 1 flavor gauge field configurations currently gen- 
erated by ETMC [2U]. Finally, one should use the obtained BB potentials 
as input for phenomenological considerations to answer e.g. the question, 
whether two B mesons are able to form a bound tetraquark state. 
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